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Abstract 

We analyse polarised Drell-Yan process using the factorisation method and 
derive operator definitions for polarised parton distribution functions. We 
demonstrate that a factorisation analogous to that in the unpolarised Drell- 
Yan case holds in this process. We study the leading order gluonic contribu- 
tion to the first moment of polarised Drell-Yan function and show that it is 
consistent with results obtained from polarised deep inelastic scattering. 
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The post Electron Muon Collaboration (EMC) IJ era has seen a revival of interest in 
the theoretical study of the polarised structure functions of the proton. The EMC exper- 
iment (polarised deep inelastic \xP scattering) gave some surprising results regarding the 
quark spin contribution to the proton spin, namely, that it is negligibly small. The most 
appealing interpretation of this result is that not only the quarks but the gluons also can 
contribute to the proton spin, to the lowest order. This is due to the fact that the product 
of strong coupling constant, a s , and the first moment of polarised gluon distribution func- 
tion, Ag(x), is of order a° @, [|. In the parton model, using the electromagnetic probe 
to evaluate the gluonic coefficient functions to leading order, one needs a prescription to 
regulate mass singularities present, and depending on the prescription the gluons may or 
may not contribute to the first moment of the polarised proton structure function, gi P (x) 
[§, || . Even if one uses the weak probe to evaluate the gluonic coefficient function, one 
would draw similar conclusions for the moments of the weak structure functions #1,3,4(0;) 
[j|]. On the other hand, in the operator product expansion || and factorisation method 
[0. §, the gluonic coefficient function is free of mass singularities and its first moment 
vanishes identically giving no gluonic contribution to the first moment of gi P (x) to or- 
der a s ||. Hence, using the factorisation method, in a class of deep inelastic scattering 
(DIS) processes probed by electromagnetic |§ and weak currents ||, one finds that the 
gluonic contribution to the first moments of polarised structure functions, g% (x) and 
#1,3,4(0;), is zero. Under these circumstances, it is therefore instructive to analyse the po- 
larised Drell-Yan (DY) processes as possible sources to measure the gluonic contribution 
to proton spin. Phenomenologically UI0 |, this process has drawn much attention after the 



EMC experiment. Further, the experimental possibilities [IJJ are also bright with the 
recent technical developments in the acceleration of polarised protons, using the so-called 
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Siberian snakes, which facilitate the polarised proton-proton collision at collider energies. 

In this letter we analyse the polarised DY process using the factorisation method and 
show that the polarised DY distribution function W(r, q 2 ) can be expressed in terms of 
known operators available in the theory. Further including the appropriate gluonic oper- 
ators we write a factorisation theorem for polarised DY process and show that the mass 
singularities factorise leaving the gluonic coefficient functions free of mass singularities to 
order a s . 

Consider the longitudinally polarised DY process Pi{piY + ^2(^2)^ — ► 7*(?) + X(Px), 
where 7* is the virtual hard photon. The cross section of the process is 

* = alf^ - ^ / A 6{q2 ~ Q2) e{q0) (1) 
where K is the flux factor, Q 2 is mass of the virtual hard photon, the sum over pho- 
ton polarisations yields the factor — g^ v and W^ v is the Fourier transform of product of 
hadronic currents sandwiched between two polarised proton states. The polarised DY 
function Wp 1 p 2 (t, q 2 ) is defined as 

W Pl p 2 (T,q 2 ) = - JL / d 4 q 5{q 2 - Q 2 ) 9(q°) J d 4 x e"** 

< p 1 s 1 ;p 2 s 2 \ J „,{x) J M (0) |pisi;p 2 s 2 >s P ■ (2) 

We choose the centre of momentum light cone frame where pi = (p±, 0, 0j_), p 2 = 
(0,p2;0_i_), q = (q + ,q~ ,0j_), x = (x + , x~, and the light cone variables are defined 
asp 1 * 1 = (p° ±p 3 )/\^2. Using the DY limit pt,p 2 — >■ 00, q + ,q~ — > 00 with r = q + q~/pfp2 
(r = Q 2 / s, where s is the cm energy) fixed, we can express the product of hadron currents 
in terms of the parton currents which are made up of parton field operators coming from 
two different protons. In order to proceed, we need to express the current-current product 
so that the dependences on the two proton momenta separate out. We do this by a Fierz 
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transformation whereby the product of local bilinear (with respect to two different proton 
states) operators is expressed as sum of bilocal linear operators. Hence the DY function 
can be neatly separated out in terms of two proton states and we obtain, 



-iq^x g— iq x n 

247T J q+ q~ ~ y ~ 2q+q~ ' 
[<PjSi| V a (x + ,x") 7 M 7 5 ^a(0) \p 1 s 1 > x 

< P2S2I Vb(0) 7 M l^b{x + ,x~) \p 2 s 2 > +(a «-> 6) , (3) 

where i runs over different quark flavours, and the scalar, vector and pseudoscalar terms 
which appear in the Fierz identity do not contribute for the polarised DY process. From 
the asymptotic behaviour of Fourier transform of singular functions, it is easy to see that 



the dominant contribution comes from the region near the singularities |12|] . Hence, using 
q + = x a pf and q~ = x b p 2 , we have 

dx n dxh 27re? t , 



■ J Xn Xh <J 



1 f 

— J dx~e~ iXaP i x ~ < p x si| ^ a (x~) 7 + 75^a(0) > x 
1 

47T 



cb + e iXbJ, 2 x < p 2 s 2 \ ip b (0) 7 75 ^b(x + )\p 2 s 2 > +(a 6) 



,(4) 



where we have also used the fact that < p ± s\ ip 7 T 75 ^ |p ± s >= 0. The above expression 
provides neat separation of short and long distance physics. By noting that the Dirac 
delta function corresponds to the q, q — > 7* subprocess, the terms in the bracket can be 
identified with the polarised quark and antiquark distribution functions of the proton, 
which are consistent with the DIS polarised distribution functions || and hence confirms 
the universality of the parton distribution functions. 

To make the definitions of the parton distribution functions gauge invariant, one has 
to insert the path ordered exponential given below. With this, one can generalise the 



above expression by replacing delta function by hard scattering coefficient function H ab 
and adding appropriate gluonic operators of the theory. We thus have 



W Pl p 2 (r,q 2 ) = I — f — fAa(x a ,fi 2 ) H ab (^—,^-) fAt(x b ,fi 2 ) + 
7fi Jr x a J ^ Xb 1 XaXb V- 2 

(a~b)+---, (5) 

where fi 2 is the factorisation scale, Wp 1 p 2 = s(W++ — W+-), /ao/a is the polarised 
probability distribution of the parton a(= qi,(ji,g) in the polarised hadron target A, 
defined as f/\a/A + = {f a+ /A+ — fa-/A+) an d H ab is the polarised hard scattering coefficient 
defined as H ab = \(H a+b+ — H a+b _) where the +(— ) sign denotes the polarisation of 
the proton or parton along (opposite to) the beam axis. The ellipses in eq.([5]) represent 
the higher twist contributions. The gauge invariant definitions of polarised quark and 
antiquark distributions in the proton are identified using eq. (f| ) in terms of the connected 
matrix elements of the bilocal operators as 



Air 
1 



dCe 



< p, s\ tf o (0, T, 0±) 7 + 75 Ql * 6 (0) \p, s > c , (6) 



f±i(x) = - I <IC < " " 



< P ,s\ * a (o) 7 + 75^ at * 6 (o,r,o±) \p,s> 



(7) 



In terms of the appropriate gluonic operators in the theory the gluonic distribution func- 
tion is defined as 



fAg(x) 



Airxp^ 



-ixl; 



< P}S \ f+»(o,c,o±) Qt K b (°) b> s >< 



<p,s\ F^(0) Qt F+*(0,r,0±) lp,a> t 



(8) 



where F£ v is the gluon field strength operator, Q a b = Vexp [ig Jq d( A + (0,( ,0±)] b , V 
denoting the path ordering of the gluon field operators A^. For known targets like quarks 



and gluons the distribution functions are normalised as 

f^+A^iz) = h 5(1 -Z) 5 a ,(q,q), (9) 

a(h) 

fM.(z) = h5(l-z)5 a , g , (10) 

a(h) 

where h = ±1 denotes the helicity of the incoming parton. In general, the above matrix 
elements cannot be calculated within the domain of perturbative QCD, owing to the com- 
plex nature of the proton target. But these are calculable for specific parton targets like 
quarks and gluons. Leading order corrections to these distributions are both ultraviolet 
(UV) and mass sensitive. These can be regularised and renormalised in any scheme. We 
use dimensional regularisation to regularise the UV divergences and MS scheme to renor- 
malise these distributions. This choice is preferred as it ensures gauge invariance and 
relativistic invariance although it destroys scale invariance. The mass singularities are 
avoided by considering on-shell quarks in the cut quark loop and off-shell gluons. Unlike 
the distribution functions f^a/A, the hard scattering coefficients are determined in the 
theory without any reference to the target. These coefficients may also be mass and UV 
sensitive and should be regularised in the same scheme. 

The n th moment of the DY function can be related to the n th moment of the distri- 
bution functions and the hard scattering coefficients as 

W AB (n,q 2 ) =£/ ¥ (n,/i 2 ) H ah (n,^) /A,(n,/i 2 ), (11) 

where f(n) = Jq x n ~ 1 f(x)dx. Here we can easily identify /ao/a^, A* 2 ) with the matrix 
elements of the twist two local operators. 

One of the interesting offshoots of the EMC interpretation is that a s Ag is of order 
Hence the leading order contribution to the DY process comes from the following sub- 
processes qq — > 7*, g(q, q) — > j*(q, q) and gg — > 7*gg. We therefore use the factorisation 

5 



theorem to calculate the gluonic contribution to the spin dependent DY function to order 
a s and the quark contribution to lowest order. To this order, the quantity of interest is 



W PlP2 {T,q 2 ) = E/ ~!r It -?fMi/Pi( x a,^ 2 ) H <im(-r— A) fAqi/Pzixb, n 2 ) + 

l Jt X a J — Xb X a Xb [I 

E / < ~~r~ I T ~~ /a»/a(z«,a* 2 ) H qig (-^—,^-) f Ag /p 2 (x b ^ 2 ) + 

I J T X a J— Jj\, X a X() /J, 

E / — f T — /Aft/Pi (X a ,// 2 ) Hq. g (-^—,^-) f A g/P 2 (x b ,fi 2 ) + 

■ Jt X n J— Xh XnXh U 



Using the Born diagram (fig.l), we get 

~ 2vre 2 

H qq -(r) = -^8(l-t). (13) 

To evaluate the gluonic contribution to the DY function Wp 1 p 2 (r, q 2 ) we replace the proton 
targets by quark and gluons and use eq.(^|), eq.(^), eq.([8]) and eq. (|12]) and express the 
resulting 'compton' subprocess to order a s as 

W qg (t,q 2 ) = I — L — Uqi/giVa, l?) H qiqi {-^—^) fA qi / q {Vb, V?) + 

i j t Va J — Vb VaVb r 1 

E / ^ L IT fMidVa,^) Hqig{--—> 3) fAg/giVb, fJ?), (14) 

where m is the mass of the quark, k 2 is off-shell mass of the gluon and f (= rjx a Xb) is the 
DY variable in the parton level. The only unknown in eq.fll4]) is f^q/g to order a s , which 
is evaluated using eq. fl7|) in the dimensional regularisation scheme. The only contribution 
to f A q/ g (h) comes from the cut triangle diagram (fig. 2) with quark and antiquark loops. 
The matrix element is only UV divergent if we keep quarks on shell and gluon off-shell. 
On regulating we get 



,/Wf,/4) = h^- 

1(h) 47T 



/4 1 , ™ 2 



^-"' V-m-fi -'^-mi-f) 



(15) 



where p? R = {j, 2 exp(jE + -) and h is the helicity of the gluon target. Since H qq is given 
in eq.(|l3|), H qg can be found if W gg is calculated. As the compton subprocess diagram 
(fig.3) has only mass singularities, we keep all the masses non zero. Neglecting the power 
corrections in the large q 2 and s limit, we get 



W qg {t,q 2 ) = £> s -g 



p 2 r « 2 m - t) 2 i 

' '2r-l)ln 2 gU 7 -^(3f 2 + 2f-3) + 

T{m z — k l r[\ — T)) 2 

2 

m 



m? — k 2 r(l — f) 

Now we can evaluate the hard scattering coefficient, by using the above results 



(16) 



~ e 2 
H qg (f,q 2 ) = -J2 a sf 



(2f - 1) (ln n T + In^f) + l(3f 2 + 2f - 3) - 1 



:i--)- "V 2" >r ™ 



= ^- 9 (r,g 2 ). (17) 

We find that a similar analysis of the weak DY processes also yields the same quark and 
gluon hard scattering coefficient functions but for weak coupling factors. Note that the 
above expression is free of any mass singularities, verifying the proposed factorisation to 
order a s in the gluonic sector. In addition, since the first moment of H qg (f,q 2 ) vanishes, 
there is no gluonic contribution to the first moment of DY function through the 'compton' 
subprocess for both electromagnetic and weak processes. 

In the parton model analysis of the above problem, the calculation of gluonic coefficient 
function in the DY function involves only W qg (f, q 2 ) term. This being mass singular, a 
particular choice of prescription has to be made to regulate these singularities. Depending 
on the prescription one gets different results: for m 2 ^ 0, k 2 = the first moment of 
W qg (f, q 2 ) is zero, while for m 2 = 0, k 2 ^ we have, 

dtW qg (r,q 2 ) = -U^f\. (18) 
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It is interesting to note that one arrives at the same conclusion even in the polarised 
DIS case. Invoking universality of parton distributions we expect the first moment of 
Wp 1 p 2 (r, q 2 ) in the parton model using the prescription m 2 = , k 2 ^ 0, to be of the form 
Aq' Aq', where Aq' = Aq — (a s /47r)Ag ||. Note that eq. (|l8|) verifies this to order a s . 
This implies that the first moment of W gg (f,q 2 ) is (Aire 2 / 3) (a s /47r) 2 . Assuming this to 
be true in the proposed factorisation scheme, one finds that the first moment of H gg (f, q 2 ) 
is also zero. Explicit computation of H gg involves evaluation of eight diagrams and three 
body phase space integral. We will be presenting these details elsewhere and expect that 
it will corroborate that the first moment vanishes. Hence, in the factorisation scheme, 
gluons can not contribute to the first moment of polarised DY function. 

We have studied the polarised DY process using in the factorisation method and 
provide operator definitions for the polarised parton distribution functions, which are 
consistent with that of DIS, hence confirms universality. We have explicitly shown the 
factorisation of mass singularities to order a s in the gluonic sector. The gluonic contri- 
bution to the first moment of the polarised DY function from the 'compton' subprocess 
turns out to be zero. Invoking universality we conclude that the gluonic contribution to 
the first moment of polarised DY function vanishes. 

We are grateful to G. Date, D. Indumathi, M.V.N. Murthy, R. Ramachandran and 
S. Umasankar for useful discussions and suggestions. Further we thank DI, MVNM, RR 
and SU for carefully reading the manuscript and making useful suggestions. 
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Figure Captions 

Fig 1. Quark Born diagram. 

Fig 2. Triangle diagrams contributing to polarised distribution functions f& q / g (h)- 
Fig 3. 'Compton' subprocess. 
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